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Abstract 

Using Dyson-Schwinger equations, the nucleon propagator is 
analyzed nonperturbatively in a field-theoretical model for the 
pion-nucleon interaction. Infinities are circumvented by using 
pion-nucleon form factors which define the physical scale. It is 
shown that the correct, finite, on-shell nucleon renormalization 
is important for the value of the mass-shift and the propagator. 
For physically acceptable forms of the pion-nucleon form factor 
the rainbow approximation together with renormalization is in- 
consistent. Going beyond the rainbow approximation, the full 
pion-nucleon vertex is modelled by its bare part plus a one-loop 
correction including an effective A. It is found that a consistent 
value for the nucleon mass-shift can be obtained as a consequence 
of a subtle interplay between wave function and vertex renormal- 
ization. Furthermore, the bare and renormalized pion-nucleon 
coupling constant are approximately equal, consistent with re- 
sults from the Cloudy Bag Model. 
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1 Introduction 



The interest in a consistent analysis of the effects of the pion cloud on nucleon 
properties is nurtured by two areas of recent research. For one, the extrapo- 
lation of nucleon lattic data, quenched or unquenched, to physical values of 
the quark (or pion) mass is mainly determined by pionic effects |]T], 0|. Fur- 
thermore, these effects will pose constraints on the development of covariant 
nucleon models such as those presented in Refs. [0, ^. The most basic ob- 
servable of interest is the nucleon mass or, more precisely, the nucleon mass 
shift associated with the pion cloud. In a previous study f^, the connection 
of covariant Euclidean nucleon mass shift calculations to extant results from 
the Cloudy Bag Model (CBM) g and Chiral Perturbation Theory (xPT) 
1^, 1^ was analyzed. A nonperturbative analysis using the Dyson-Schwinger 
(DS) equation in rainbow approximation (with no renormalization) suggested 
that almost all of the nucleon mass shift can be attributed to the covariant 
one-loop pion dressing. We will extend this analysis by supplementing the 
DS equation with correct on-shell renormalization conditions (Section 
To properly implement those, one has to forgo an angular approximation 
employed in Ref . |^ . In Section |^ we calculate a one-loop correction to the 
ttNN vertex which is subsequently employed in the DS equation. In the last 
section, we summarize and present our conclusions. 



2 The model 



We consider a pseudovector Lagrangian for the vrA^ interaction which reads 
in Euclidean space 

C^NN = 27^75 TVl> . (9^7r) . (1) 

z,IVl]\f 

The DS equation for the nucleon propagator G is given by 

G-\p) = G^\p)+ii>T.y{p^) + T.s{p^), (2) 

ii>A{p^)+B{p^) = Z^iii^ + ZmM) + Z, j ^T"^ G{k)T D{p-k) , 

r° = {i^-}^)i,Tg. (3) 

Because of the small mass of the pion, we can approximate the full pion 
propagator D{q) by the free scalar propagator, D^{q) = (g^ + mD~^. r° and 
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r stand for the free and the full ttNN vertex. The renormalization constants 
Zi, Z2 and Zm refer to the nNN vertex, nucleon wave function and nucleon 
mass respectively. The relation to bare quantities is given by 

Z2 Mbare , 
9 = ^ ^bare , Mn = ^ • (4) 

To account for the compositeness of the particles, we introduce a form factor 
at each nNN vertex into this idealized field-theoretic model, 

rM(p,A;) ^ T^^\p,k)u{{p-ky) . (5) 

In the case where both nucleons are on shell (p^ = k"^ = — M^), this form 
factor is most naturally related to the axial form factor of the nucleon. It is 
commonly parameterized as a dipole, 

A — ^ 

but we will also investigate (as in Ref. 0) the exponential form 



^ir) = (T^-f , (6) 



2n f + ml 



= exp y--^ ) ■ (7) 

From neutrino scattering experiments, the dipole width parameter (in the 
on-shell case) is determined as A = 1.03 ± 0.04 GeV [jlO|- This sets roughly 
the scale for our considerations. 



Even after the introduction of form factors, we assume that it is sensible to 
extract nonperturbative information from the DS equation (^. We note that 
now the loop integral in the equation is convergent and all renormalization 
constants will therefore be finite. A commonly used approximation to it is 
the rainbow approximation F = F'^. The full irNN vertex F fulfills its own 
DS equation which takes the symbolic form: 

F = Zi F° + F° (GWgI'I) K4 . (8) 

Here, K4 is the full, off-shell, amputated N — N scattering matrix. One sees 
that the rainbow approximation neglects the last term and enforces Zi = 1. 



We employ on-shell renormalization for the unknown nucleon propagator, 
G. This requires that G has a pole at = — (the physical mass), 

M^(Z2 + Sv(-M^)) = S5(-M^) + ZsZmM^ , (9) 

Mbare-M^ = M^(ZA,-1) = ^ (M^Sy (-M^,) - (-M^,) ) . ( 10) 



3 



The renormalization constant Z2 is determined by the condition that the 
residue at the pole be unity, 



-Ml) 



V 



Qp2 



-2M 



N 



1 ^ (11) 
Z2-I. (12) 



Mass renormahzation ensures that the self-energy has cuts starting from the 
physical thresholds. Furthermore one can show that the spectral densities 
of the solution G multiplied by Z2 are properly normalized (i.e. Z2 can be 
interpreted as the probability of finding a "bare" nucleon inside the pion- 
dressed one). 



In the exploratory study |^ the DS equation, Eq. (^, was solved in rainbow 
approximation after putting Zi = Z2 = Zm = 1 and disregarding the above 
renormalization conditions. The mass of the dressed nucleon was then found 
as the solution of -M|, A(-M|,) + 5(-M|,) = 0, < M^. In this 
simplified scenario it was possible to employ a certain angular approximation 
to the loop integral in eq. (0) to calculate all angular integrals analytically. 
However, if one incorporates the renormalization conditions (|10]) and (|12D, 
this angular approximation fails because it underestimates the slopes of Sy 
and Tjs considerably (and these enter the expression for Z2). As a result, we 
must resort to a numerical computation of one angular integral. Since we 
have to evaluate the renormalization conditions on the nucleon mass shell, it 
is necessary to continue the DS equation to complex momenta. This intricate 
procedure is outlined in Appendix ^ 

We illustrate the results for the mass shift in Fig. |^, employing (as in Ref. 
0) as the renormalized coupling constant g = Mjq / f-K, i.e. qa = 1- We chose 
a dipole form factor with cut-offs A in a range compatible with the measured 
axial form factor. [| Indeed, the mass shift for the unrenormalized rainbow 
treatment and the one in one-loop approximation are very close. In contrast, 
the mass shift in the renormalized rainbow treatment is larger by a factor 
which rises from 1.8 (A = 0.96 GeV) to 2.9 (A = 1.2 GeV). The explanation 
for this peculiar behaviour can be found in the right panel of Fig. |I| where the 
ratio of bare to renormalized coupling is plotted. This ratio equals I/Z2 in 
the rainbow treatment and is considerably larger than 1. If one divides the 
DS equation, Eq. (||), by Z2, one sees that the loop integral is proportional 

"'^For technical reasons we consider only A > Afjy. Otherwise poles from the form factor 
enter the integration domain, adding technical complications which are of no interest here. 
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Figure 1: Left panel: nucleon mass shift in the renormahzed treatment (de- 
fined as 5M = Mtv — Mbare) and in the unrenormahzed treatment (defined 
as SM = Md — Mn) of the rainbow approximation. Right paneh the ratio 
gh&rc/ g = Zi/Z2. 



to the bare couphng and therefore drives the mass shift to larger absolute 
values. 

A repetition of the calculations with the physical qa = 1-26 only aggravates 
the difference between renormalized and unrenormahzed results. Whereas 
the mass shift in the unrenormahzed rainbow treatment scales with g\ (as 
the one-loop result), the mass shift in the renormalized case shoots up to 
values between 400 MeV and 1.9 GeV (depending on A), because of the 
nonlinear nature of the DS equation. 

The question arises whether such a strong renormalization of the pion- 
nucleon coupling constant, as visible in the rainbow solutions, is physically 
reasonable. An analysis of pionic corrections to nucleons in the CBM to 
one-loop order P, |ll| reveals that Zi ^ Z2 (and therefore (?barc ~ g) if 
one includes the A resonance in the analysis of the self-energies and the 
pion-nucleon vertex.^ This is a strong indication that the rainbow approx- 
imation is unreliable for this problem: there Zi = 1 is set artificially and 
Z2 ~ 0.6 . . . 0.7 for physical irNN form factors. 

Therefore, we will model the full ttNN vertex F by the bare one plus a one- 



^Thc analysis there corresponds to a non-relativistic expansion of positive-energy 
graphs in time-ordered perturbation theory. 
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loop correction that includes the A in an effective manner, i.e. as a spin-3/2 
particle described by Rarita-Schwinger spinors. Ideally, we would like to 
employ self-consistent propagators of and A in this study, but solving a 
coupled system of two-loop DS equations for and A is beyond the scope of 
this study. Instead, we pursue a modest but nevertheless resource-consuming 
modification: the one-loop vertex correction to F will be calculated using free 
and A propagators and the result will be inserted back into the DS eq. @. 
By that token, the problem becomes numerically tractable, and we expect 
the bulk of the effect on the solutions to be buried in a reasonable estimate 
for Zx. 



3 Corrections to the ttNN vertex 




Figure 2: One-loop corrections to the ttNN vertex. 



We calculate the covariant one-loop correction to the irNN vertex by includ- 
ing the A as an effective degree of freedom. Pictorially the equation for the 
vertex is displayed in Fig. ^, and as mentioned above we take for the nucleon 
and A propagators free spin-1/2 and spin-3/2 propagators, respectively, 

^^P) - p2 + ' ^^^^ 

3 ^ ^ 3 Mi 3 Ma ^ ^ 

with Mjv = 0.94 GeV and Ma = 1.23 GeV are the physical masses of both 
particles. For the vrA^A and vrAA interactions we employ tree level vertices 
derived from the simplest covariant interaction Lagrangians, i.e., 

/:.AA = f^^^^7'^75T^;^v,..(9,.)^^0AA = f^^75T^;^(17) 
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Here q is the momentum of the incoming pion. For the isospin 1/2-3/2 
transition matrix we adopt the convention 



^3/2 _ 

1/2 ~ ^\m,lk 



As indices in the Clebsch-Gordan coefficient, M, m, k stand for the isospin-z 

-,3/2 

3/2 



Q In 

components of A, and tt. Furthermore, the isospin matrices TJ,^ are just 



the ones for the 4-dimensional SU (2) representation. 

For the numerical values of the coupling constants ^fyrAfA and ^fyrAA we relate 
them to g via SU{Q) quark model expressions. This of course leaves room to 
fine-tuning which is not the main interest here. To apply the quark model, 
we note the non-relativistic limit of both vertices, 

pO ^ 2z A'ttAA ^3/2 , jj3/2 X ,„„x 

i.AA ^ ^^3/2 1-^3/2 • i^Uj 

The transition matrix S^yg identical to T'^i^^-, it just refers to the spin degree 
of freedom. Comparing to the expressions in Ref . , we readily find 

/72 

S'ttTva = 9 , (21) 



g^AA = - g ■ (22) 



Here, we employ for the renormalized coupling constant the physical value 
g = gAMN/fiT with gA = 1.26. Having fixed the strength of our interactions, 
we proceed now to the calculation of the renormalization constant Zi. We 
choose to fix it at the virtual point where both nucleons and the pion are on 
shell, i.e. 

m^uip) = ^u{k){l^-'i>)^,u{p) (23) 

[p2 = ^2 ^ _ ^ _^2] ^24) 

This has the advantage that the uNN form factor should reduce to 1 at this 
point. 

The resulting Zi as a function of the dipole cut-off is plotted in Fig. |. From 
the right panel one can see that the loop-nucleons contribute less than 0.05 to 
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Figure 3: Left panel: Vertex renormalization constant Zi. The different 
contributions to 1 — Zi, stemming from the loop diagrams shown in Fig. 
are labelled with Vmn with m^n ^ {A^, A} labelling the baryons in the loop. 
Right panel: Zi, Z2 and their ratio in one-loop approximation. 



\ — Z\. If this were the whole story, the rainbow treatment of the DS equation 
would seem to be justified, yielding mass-shifts > 500 MeV as demonstrated 
before! The bulk of the difference \ — Z\ comes from the two graphs with 
one intermediate A. The resulting bare coupling is somewhat lower than the 
renormalized one. If the A were also included in the nucleon self-energy, we 
would expect a somewhat larger bare coupling since its contribution lowers 
Z2 additionally. Keeping that in mind, we expect the present results for the 
mass shift in the improved DS equation to be a lower bound. 

Our one-loop improved model for the vertex can now be inserted into the 
DS eq. (0). Technically, we proceed by projecting the Dirac structure of the 
vertex onto basic covariant matrices, 

rDirac(p,A;,/) = (Z75) \/l(p^A;^p■/c) + (J^^75) ^2(p',/c',p-fc) + 

(^T75) V^3(/,A;^p■A;) + (^^TJ^^75) Vi{p\k\p-k) ,(25) 

= k-p^, (26) 

from which the scalar functions Vi,,a can be readily traced out, using the 
program FORM |]T2|. The three remaining scalar integrals are evaluated 



numerically. We calculate the functions Vi on a three-dimensional grid in 
the variables p^, fc^ and p ■ k, needed for the solution of the DS equation. 
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Figure 4: Left panel: Mass shift for the vertex-improved DS solution in 
comparison to the one-loop approximation. Right panel: Zi, and their 
ratio in the DS solution. 



The results for mass-shift and the wave-function renormalization constant 
Zi are depicted in Fig. ^. In the physically interesting region A G [0.95, 1.05] 
GeV the mass shift stays rather fiat, at a value around 200 MeV. For a harder 
form factor it actually drops, which can again be explained by looking at the 
ratio Qharc/ 9- beyond A ~ 1.1 GeV gharc/ 9 is becomes less than 1/2. As 
explained before, this ratio enters the pion loop integral in the DS equation 
through Zi. Since Zi itself is very low there, the one-loop treatment of the 
ttNN vertex can be questioned. 

The feature of a plateau in the mass shift for a certain range of cut-off 
parameters (which happens to coincide with the experimental results for the 
axial form factor) is actually a desired one in effective models with cut- 
off functions, since it indicates the relative independence of the results on 
the specific choice of cut-off. Whether that still holds after self-consistent 
inclusion of the A resonance is an open question. Certainly, by the arguments 
given above, its proper inclusion should lead to ghnre/ 9 > 1 and a larger 
mass shift. Keeping this in mind, a value —6M ^ — (5Mi_ioop ~ 200 MeV 
constitutes a lower bound. 

In Fig. ^ we show the results for an exponential nNN form factor. Since the 
exponential is an entire function, we can investigate also somewhat softer 
form factors without encountering analytical problems. One-loop mass shifts 
are the same for dipole and exponential, if A ~ 1.2A. The result for the self- 
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gt,^Jg, improved DS solution 

NuCleOn mass Sllift exponential form factor 




0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 
A [GeV] A [GeV] 



Figure 5: Mass shift and renormalization constants arising from a pion- 
nucleon form factor of exponential type. 



consistent mass-shift looks overall very similar to the previous case. Its curve 
flattens around A = 0.8 GeV at a value of approximately 190 MeV and then 
drops because of the rapidly decreasing bare coupling constant. 



4 Summary and Conclusions 



We have investigated the covariant Dyson-Schwinger equation for the nucleon 
in a field-theoretic model for the pion-nucleon interaction with a pseudovec- 
tor Lagrangian. Its treatment in the commonly used rainbow approximation, 
including proper renormalization, leads to a strong renormalization of the 
pion-nucleon coupling constant, gharc/ 9 > 1, in contradiction to perturba- 
tive one-loop calculations. We were therefore led to calculate the one-loop 
perturbative correction to the vertex, including the A resonance, and to re- 
solve the DS equation. For physically reasonable values for the cut-off in 
the dipole pion-nucleon form factor, A ~ 1 GeV, we find 5'barc/fi' ~ 1 and a 
rather stable value for the nucleon mass-shift of around —200 MeV, consis- 
tent with one-loop results in both covariant treatments and semi-relativistic 
approaches such as the cloudy bag model. The fully self-consistent inclusion 
of the A resonance is expected to raise the absolute value of the mass shift 
even further. 
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A Solving the Dyson— Schwinger equation for 
complex momenta 

Here we present a more detailed account of the procedure to solve the DS 
eq. (PD, 



ii>A{p')+B{p') = Z2{t^ + ZMM)+t:^Ev{p') + ^sip') . (27) 



z^Ey{p') + J:s{p') = Z,I^V\p-k)x 



^ T{p,k) D^{p-k) . (28) 



ilj: A{k'^) + B{k 

In order to evaluate the renormalization conditions (p^p!2D , we will need (as 
shown below) the functions AijP') and BijP') for p^ G [—{Mm — rn.„y, oo) - 
i.e. if p = (0,^4), p4 may assume imaginary values. 

We introduce the two component vectors S(p^) = (Ss'(p^), Sy(p^)) and 
S'(p^) = {A{p'^) , B [p'^)) . After applying suitable traces and doing the two 
trivial angular integrals we arrive at an equation of the form 

3 f°° k'^ dk'^ 

^'^^'^ ^ (2^0^ Jo k^A^mTB^k^) ^^^^ 

\/l — z"^ dz 



I ,^,2! ft IM K,Ap\k\z)S,{k' 
J -I p^ + k^ + — 2p^\k\z 



where z = cosip refers to the angle between the Euclidean vectors p and k. 
If we want to evaluate the loop integral for < (p4 imaginary), we have 
to note that for p^ + < poles from the pion propagator cross the real k^ 
axis. These have to be avoided by a suitable deformation of the integration 
contour. This situation is depicted in Fig. p. Alternatively, the original real 
path may be retained provided that a loop contour around the poles is added. 
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integration contour 
Re 



Figure 6: The ^4 integration in the pion loop. Once the singularities have 
crossed the real axis, the proper integration contour consist of the real axis 
and a loop contour around the singularities that have crossed the axis. 



This leads to 



(27r)3 
1 



a/1 — 2^ dz 



X 



1 + /c^ + — 2^4/^42; 



Ki,{p\k\z) S,{k') + 



K^,{p\k\z) 
k^A^{k^) + B^{k^) 



poles 2 \/ 
S.ik') 



X 



(30) 



(p— fc)2=— ml 



The loop contour in k^ picks up residues (— 27ri) / (— 2i k^ + m^) from the 
pion propagator, and the integrand for the remaining integral in 3-space has 
to be evaluated at the pion poles (thus the value of z to be used in Kij is 
determined). We convert the integral over k into an integral over the squared 
four-momentum k"^ using 



dk' 



\k\ 



de k'' F(A:^p2^ 



2y/k' + ml 



F{k\p' 



(31) 
(32) 



Hence, we arrive at the final expression for the self-energies, valid also for 
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timelike momenta up to —p^ = M^: 



3 r dk^ 



(27r)3 7o PA2(A;2) + 52(F) 

J_l p2 + fc^ + 777,2 - 2P4/C42; 



K,,{p\k\z) 
k^A\k^) + B^{k^] 



Sjik^ 



(33) 

(p— fc)^=— 



Now one sees clearly that in order to evaluate the renormalization conditions 
at p2 = — M^, one needs to know the self-consistent solution to v4(p2) and 
i?(p2) for the interval p2 g [_(j^f^ _ 777^^)2^00). The DS equation can now 
be solved iteratively, performing both z and k'^ integrations numerically. By 
virtue of the pion-nucleon form factors, the numerical treatment is not ham- 
pered by ultraviolet divergences and stable results are achieved by using a 
mininum of 100 mesh points for each integral. 
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